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Abstract 

Energy dependence of space-time parameters of pion emission region at 
freeze-out is studied for collisions of various ions and for all experimentally 
available energies. The using of femtoscopic radii scaled on the averaged 
radius of colliding ions is suggested. This approach allow/s the expansion of 
the set of interaction types, in particular, on collisions of non-symmetrical 
ion beams w/hich can be studied w/ithin the framew/ork of common treatment. 
There is no sharp changing of femtoscopic parameter values with increasing 
of initial energy. Analytic functions suggested for smooth approximations of 
energy dependence of femtoscopic parameters demonstrate reasonable agree¬ 
ment with most of experimental data at ^/snn > 5 GeV. Estimations of some 
observables are obtained for energies of the LHC and FCC project. 


1 Introduction 

At present femtoscopic measurements in particular that based on Bose-Einstein 
correlations are unique experimental method for the determination of sizes and 
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lifetimes of sources in high energy and nuclear physics. The study of nucleus-nucleus 
{AA) collisions in wide energy domain by correlation femtoscopy seems important for 
better understanding both the equation of state (EOS) of strongly interacting matter 
and general dynamic features of soft processes. The discussion below is focused on 
specific case of femtoscopy, namely, on correlations in pairs of identical charged pions 
with small relative momenta - HBT-interferometry - in nucleus-nucleus collisions. 
The general view for phenomenological parameterization of correlation function (CF) 
for two identical particles is discussed in the [HE]- Below the experimental results 
obtained for AA collisions within the standard 3d approach are taken into account 
[2]. The set of main femtoscopic observables Qi = = {X, Rs, Ro, R\} is 

under consideration as well as the set of some important additional observables 
which can be calculated with help of HBT radii Q 2 = {^^ 2 }^=! = {Ro/Rs,S,V}. 
Here 6 = Rl-R^V = {2Trf^^R‘^Ri is the volume of source at freeze-out. The set 
of parameters Qi characterizes the correlation strength and source’s 4-dimensional 
geometry at freeze-out stage completely. The most central collisions are usually 
used for study the space-time characteristics of final-state matter, in particular, 
for discussion of global energy dependence of femtoscopic observables. Therefore 
scaled parameters QI, i = 2 — 4, S and are calculated as follows mm- 

= RJRa, i = s.o.l; <5" = 5/Rl-, F" = V/V^. (1) 

Here i?A = RqA^^^= 4ttR\/‘4 is radius and volume of spherically-symmetric 
nucleus, i?o = (1-25 ± 0.05) fm [21|4]. The change i?A —t {R^) = 0.5(i?Ai + Ra^) 
is made in the relation (Q in the case of non-symmetric nuclear collisions [HE]. 
In general case the scale factor in 0 should takes into account the centrality of 
nucleus-nucleus collisions. The normalization procedure suggested in [J allows the 
consideration of all available data for nucleus-nucleus collisions |2]. As development 
of previous analyses [HE] the proton-proton [pp) results at high energies [5||^ are 
also considered here with replacing i?A —> in (Q. 

2 Energy dependence of space-time extent of 
pion source 

Detail study for (quasi)symmetric heavy ion collisions [HE] demonstrates that the 
fit function (e = snn/sq, sq = 1 GeV^) 

fWsNN) = ai [1 -F a2(ln£)“3] (2) 

agrees reasonably with experimental dependence 5](-\/s)vw), f = 1 — 4 at any 
collision energy for A and at \/snn > 5 GeV for HBT radii. Fig.shows the en¬ 
ergy dependence of A (a), scaled HBT-radii (b - d) and Ro/Rs ratio (e) for both 
the symmetric and non-symmetric collisions of various nuclei. Fits of experimen¬ 
tal dependencies for AA interactions are made by (|^ in the same energy domains 
as well as for (quasi)symmetric heavy ion collisions. The numerical values of fit 
parameters are presented in Table[2 fit curves are shown in Fig.[^ by solid lines 
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Figure 1: Energy dependence of A parameter (a), scaled HBT-radii (b - d) 
and ratio Ro/Rs (e) in various collisions. Experimental data are from [51[51|S]. 
Statistical errors are shown (for NA44 - total uncertainties). The solid lines 
(a - d) correspond to the fits by function ([^ and dashed lines - to the fits by 
specific case of (H) at fixed 03 = 1.0. Smooth solid and dashed curves at (e) 
correspond to the ratio Ro/Rs calculated from the fit results for i?" and ii" in 
AA, dotted line is the level Ro/Rs = 1. 


for Q and by dashed lines for specific case of fit function at 03 = 1.0 with tak¬ 
ing into account statistical errors. There is dramatic growth of y^/n.d.f. values 
for fits of A data (Fig.[^) despite of qualitative agreement between smooth ap¬ 
proximations and experimental A values for range 10 < \/s//~n < 200 GeV. The 
fit by ([§ underestimates the A value at the LHC energy ^/sNN = 2.76 TeV sig¬ 
nificantly. The A values for asymmetric nucleus-nucleus collisions at intermediate 
energies \/snn 20 GeV agree well with values of A in symmetric heavy ion col¬ 
lisions at close energies. On the other hand the A for Cu-I-Cu collisions is smaller 
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Table 1: Values of fit parameters for A A data with statistical errors 


HBT 

parameter 


Fit parameter 


ai 

02 

03 

xVn.d.f. 

A 

1.21 ±0.09 

-0.30 ±0.04 

0.38 ±0.04 

3656/29 


0.717 ±0.003 

-0.051 ±0.001 

1.0 (fixed) 

3786/23 

on 

Jls 

0.656 ±0.002 

(6 ±3) X 10-5 

3.11 ±0.19 

195/25 


0.599 ±0.003 

0.019 ±0.001 

1.0 (fixed) 

280/26 

zpn 

Jto 

0.10 ±0.02 

6.3 ± 1.7 

0.068 ±0.006 

402/25 


0.758 ±0.004 

0.008 ±0.001 

1.0 (fixed) 

415/26 

Rl^ 

0.022 ±0.002 

23 ±3 

0.258 ±0.005 

502/25 


0.634 ± 0.004 

0.043 ± 0.001 

1.0 (fixed) 

615/26 


systematically than A in Au+Au collisions in energy range ^/snn = 62 — 200 GeV 
(Fig.[^). New experimental data are important for verification of the suggestion 
of separate dependencies X{^ysNN) for moderate and heavy ion collisions. Also the 
development of some approach is required in order to account for type of colliding 
beams in the case of A parameter and improve quality of approximation. Smooth 
curves for normalized HBT radii and ratio Ro/Rs are in reasonable agreement with 
experimental dependencies in fitted domain of collision energies ^/snn > 5 GeV 
(Figs.[^ - e). Dramatic improvement of the fit qualities for scaled HBT radii at 
transition from the data sample with statistical errors to the data sample with total 
errors is dominated mostly by the uncertainty in rg leads to additional errors due 
to scaling ©■ The scaled HBT-radii in pp are larger significantly than those in 
AA collisions at close energies. Because feature of Regge theory [7] the following 
relation is suggested to take into account the expanding of proton with energy: 
Rp = ’"o(l + ky/aip Ine), where rg = (0.877 ± 0.005) fm is the proton’s charge 
radius |8], parameter a'-p oc Ine because of diffraction cone shrinkage speeds up 
with collision energy in elastic pp scattering [9|. The k is defined from the bound¬ 
ary condition Rp —>■ l/m^ at e —> oo with choice of appropriate asymptotic energy 
The detail study demonstrates that the increasing of from 6 PeV 

[10] to 10^ PeV influences weakly on i?", i = s,o, I in pp collisions and calculations 
are made for the first case. The normalized transverse radii agree in both the pp 
and the AA collisions (Figs.[^, c) at = 200 GeV with excess of i?" in pp 

with respect to the AA in TeV-region. The i?" in pp is larger than that for AA 
in domain \/snn > 200 GeV. It should be stressed that the additional study is 
important, at least, for choice of i?p(e). 

The corresponding dependencies for and V” are demonstrated in Fig.and 
Fig.[^ respectively. As well as in [T|[2 results for 7r+7r+ pairs are shown in Figs.[^- 
[^also because femtoscopy parameters from the set Qi depend on sign of electrical 
charge of secondary pions weakly. The relation Ro < Rs is observed for « 11% 
of points in Fig.|^ Detail discussion for points with <5 < 0 is in the [2]. The 
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Figure 2: Dependence S'^{^/snn) in various collisions. Experimental data are 
from [a Eli. Error bars are only statistical (for NA44 - total uncertainties). 
Dotted line is the level i5" = 0. Smooth curves are derived from 0 and the fit 
results for i?”, i?" in A A. The solid line corresponds to the hts of normalized 
HBT radii by function ([^ and dashed line - to the fits by specific case i?" oc In e, 
i = s,o. 


dependence is almost flat within large error bars in all energy domain 

under consideration. Taking into account the STAR high-statistics results m only 
one can see the indication on change of behavior of (5”(-^SArAr) inside the range of 
collision energy \/snn = 11.5 — 19.6 GeV. This observation is in agreement with 
features of behavior of emission duration (At) dependence on \/snn discussed in 
[2]. The estimation of energy range agrees well with results of several studies in 
the framework of the phase-1 of the beam energy scan (BES) program at RHIC 
which indicate on the transition from dominance of quark-gluon degrees of freedom 
to hadronic matter at ^^/snn ^ 19.6 GeV. But future precise measurements are 
crucially important for extraction of more definite physics conclusions. Smooth 
solid and dashed curves shown in Fig.j^are calculated for 5” from the fit results for 
i?” and i?Q (TableQ. The calculation based on the fit function ([^ at free 03 agrees 
reasonably with experimental points at ^/snn < 200 GeV but underestimates <5” 
in TeV-region significantly. The large errors in Fig.for strongly asymmetric AA 
collisions is dominated by large difference of radii of colliding moderate and heavy 
nuclei and corresponding large uncertainty for (Aa)- Smooth solid and dashed 
curves shown in Fi g.[3] are calculated for F" from it’s definition (Q and the fit results 
for i?", (TablellJ). The fit results for normalized HBT radii obtained with general 
function ([^ lead to very good agreement between smooth curve and experimental 
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data in TeV-region in contrast with the curve obtained from corresponding fit results 
for ([^ at 03 = 1.0. There is significant difference between pp and AA collisions 
for (5" in TeV-region (Fig.[^ and for V" at \Jsnn > 200 GeV (Fig.|^. 



Figure 3: Energy dependence of F” in various collisions. Experimental data are 
from EIEIE]- Error bars are only statistical (for NA44 - total uncertainties). 
Smooth curves are derived from Q and the fit results for i?", i?p in AA. The 
solid line corresponds to the fits of normalized HBT radii by function and 
dashed line - to the fits by specific case oc Ine, z = s, 1 . 

Estimations for A, Ro/Rs, and normalized femtoscopic parameters at the LHC 
and the FCC energies are shown in Table[^ for fits of various AA collisions with 
inclusion of statistical errors, the second line for each collision energy corresponds 
to the using of the specific case of ([^ at 03 = 1.0. All the smooth approximations 
discussed above predict amplification of coherent pion emission with significant 
decreasing of A. Uncertainties are large for estimations obtained on the basis of 
results of fits by function (|^ at free 03 . Thus values of femtoscopic observables in 
Tablej^are equal within errors for general and specific case of ([^ at ^/sNN = 5.52 
TeV (LHC) and = 39.0 TeV (FCC). 

The energy dependencies for sets Qm, m = 1, 2 of femtoscopic parameters with 
taking into account the scaling relation (Q demonstrate the reasonable agreement 
between values of parameters obtained for interactions of various ions (Figs. 0 - 11 - 
The observation confirms the suggestion [I] that normalized femtoscopic parameters 
allow us to unite the study both the symmetric and the asymmetric AA collisions 
within the framework of united approach. This qualitative suggestion is confirmed 
indirectly by recent study of two-pion correlations in the collisions of the lightest 
nucleus (d) with heavy ion (Au) at RHIC. Estimations of space-time extent of the 
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Table 2: Estimations for observables based on fit results 


y/SNN, 


HBT parameter for AA 


TeV 

A 

TDn 

Jig 

TDn 

-Tto 

R^ 

5.52 

0.16 + 0.19 

0.9+ 0.2 

0.8+ 0.3 

1.06 + 0.16 


0.091 ± 0.004 

0.792 + 0.009 

0.860 + 0.010 

1.099 + 0.013 

39.0 

0.07 + 0.21 

1.2+ 0.4 

0.9+ 0.3 

1.11 + 0.16 


- 

0.836 + 0.011 

0.883 + 0.012 

1.205 + 0.015 


Ro/ Rs 

(5" 

yn 


5.52 

0.9+ 0.4 

-0.2+ 0.6 

3.5 + 1.6 



1.086 + 0.018 

0.11 + 0.02 

2.59 + 0.07 


39.0 

0.7+ 0.3 

-0.7+ 1.1 

6 + 4 



1.06 + 0.02 

0.08 + 0.03 

3.17 + 0.09 



pion emission source in d+Au collisions at top RHIC energy jT2] in dependence 
on kinematic observables show similar patterns with corresponding dependencies 
in Au+Au collisions and indicate on similarity in expansion dynamics in collisions 
of various systems (d+Au and Au+Au at RHIC, p+Pb and Pb+Pb at LHC). The 
scaling results for some radii indicate that hydrodynamic-like collective expansion 
is driven by final-state rescattering effects |12| . On the other hand the normalized 
femtoscopic parameters allow us to get the common kinematic dependencies only 
without any additional information about possible general dynamic features in dif¬ 
ferent collisions. Thus the hypothesis discussed above is qualitative only. The future 
quantitative theoretical and phenomenological studies are essential for verification 
of general features of soft stage dynamics for different collisions at high energies. 


3 Summary 

The main results of present study are the following. 

Energy dependence is investigated for range of all experimentally available initial 
energies and for estimations of the main femtoscopic parameters from set the 
(A and radii) derived in the framework of Gauss approach as well as for the set of 
important additional observables Q 2 contains ratio of transverse radii, S and HBT 
volume. There is no dramatic change of femtoscopic parameter values in AA with 
increasing of in domain of collision energies \/snn > 5 GeV. The energy 

dependence is almost flat for the in nucleus-nucleus collisions within large error 
bars. The indication on possible curve knee at ^/snn ^ 10 — 20 GeV obtained 
in the STAR high-statistics data agree with other results in the framework of the 
phase-1 of the BES program at RHIC. But additional precise measurements are 
crucially important at various y/s^N in order to confirm this feature in energy 
dependence of additional femtoscopic parameters. The normalized some HBT radii 
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and source volume in pp are larger significantly than those in AA collisions especially 
in TeV-region. The fit curves demonstrate qualitative agreement with experimental 
AA data for A at all available collision energies and for normalized HBT radii in 
energy domain y/sj^N > 5 GeV. Smooth curves calculated for energy dependence 
of parameters from the set Q 2 agree reasonably with corresponding experimental 
AA data in the most cases. Estimations of femtoscopic observables are obtained 
on the basis of the fit results for energies of the LHC and the FCC project. For 
multi-TeV energy domain the emission region of pions in nucleus-nucleus collisions 
will be characterized by decreased correlation strength, linear sizes about 8.5 — 9.5 
fm in longitudinal direction and 7 — 8 fm in transverse plane, volume of about lO"* 
fm^. 
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